Ground-States of Two Directed Polymers 
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^ ■ Abstract 



Joint ground states of two directed polymers in a random medium are investi- 
gated. Using exact min-cost flow optimization the true two- line ground-state 



\ is compared with the single line ground state plus its first excited state. It is 

O \ found that these two-line configurations are (for almost all disorder configura- 

tions) distinct implying that the true two-line ground-state is non-separable, 
even with 'worst-possible' initial conditions. The effective interaction energy 
between the two lines scales with the system size with the scaling exponents 
\ 0.39 and 0.21 in 2D and 3D, respectively. 
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I. INTRODUCTION 



The physics of disordered systems has attracted a lot of attention due to the discovery 
that the free energy of extended objects - hnes, surfaces and so on - has singular corrections 
because of the domination of zero-temperature or ground-state effects [|l|. The paradigm 
of such systems is a directed polymer in a random medium (DPRM). In this particular 
example the object minimizes its energy which is determined by two competing forces: 
the elastic energy cost of wandering on one hand and the energy gain using energetically 
favorable pins in the environment on the other hand. The result is super-diffusive behavior, 
and constrained energy fluctuations. The phase space of the DPRM problem is very rich 
depending on the nature of the correlations in the disorder and the dimensionality. In low 
enough dimensions the physics is (at arbitrary temperatures) governed by the so-called zero- 
temperature fixed point if the noise has weak enough correlations including the uncorrelated 
case. The case with one transverse dimension becomes exactly solvable in terms of the 
roughness and energy fluctuation exponents, due to a mapping to the Kardar-Parisi- Zhang 
equation p|. The values are C2 = 2/3 and 62 = 1/3, which fulfill the exponent relation 
2C^d — ^ = 6d- In the 3(=2+l)-dimensional case the roughness exponent is approximately 
0.62. 

In this paper we study the problem of two (not necessarily directed) polymers in a 
joint random medium (TPRM) with mutual interactions and focus on the repulsive 
strong coupling limit, i.e. hard core interaction. The work is related both to the question 
of the physics of fiux-lines in high-Tc-semi conductors in the low field limit, and to the field- 
theoretical issues due to the importance of the DP interaction energy. The physics of the 
problem is in general very similar to that of the one-line case |]9|JT0|J^ but shows interesting 
twists if one tries to understand the problem in the light of individual, independent objects. 
In particular, we are going to consider by numerical, exact min-cost flow optimization com- 
putations the difference in energy between the TPRM problem, the single-line ground-state 
and the 'first excited state'. The last one is given by adhering to a hierarchical picture. 



2 



in which the first polymer is first optimized given a disorder configuration, and then the 
next one is added by applying a hard-core repulsion to the bonds already taken up. The 
procedure gives us two energies to compare with the true TPRM ground state energy E2: 
the single line ground-state energy doubled, 2Ei, and the sum of the ground-state energy Ei 
and the energy of the first excited state E[ in the single line problem. The two energy dif- 
ferences, E2 — 2Ei and E2 — Ei — E[ define an interaction energy of the two polymers. In an 
earlier paper Tang studied the TPRM in hierarchical lattices and in two dimensions with 
binary disorder. His main conclusion was, for the physically more relevant real-space case, 
that the probability for an interaction energy exactly equal to zero (with binary disorder) 
decayed much faster than expected, the exponent being -2/3 instead of the -1/3 expected 
based on single-DP geometric arguments. We study both the interaction energies discussed 
above. We also comment on the topology of the TPRM ground-state. One of the main 
conclusions of our paper is that the TPRM ground-state is non- separable at least in the 
particular geometry we use. This means that the optimization of the TPRM ground-state 
can not be done in two quasi-independent steps. 

The structure of the paper is as follows. In section two we formulate the problem and 
outline the relevant scaling exponents to be studied later. Section three discusses the nu- 
merical method. In section four we give the numerical data concerning the scaling behavior. 
Finally section five finishes the paper with conclusions. 

II. TWO DIRECTED POLYMERS IN A RANDOM MEDIUM 

The continuum Hamiltonian for the TPRM problem is written in all generality as 

Jo 

The Hamiltonian describes the physics of two elastic lines (subscripts 1 and 2) in the pres- 
ence of the random potential Vr which is sample-to-sample the same for both lines. The 
longitudinal coordinate is labeled with x while the transverse coordinate (which can be a 
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vector) is hi or /i2- In the following we shall consider only two 'identical' lines, that is the 
line stiffnesses Fj are taken to be finite and equal. The random potential Vr describes point 
disorder and therefore the correlator {Vr{x, h)Vr{x', h')) oc 5{x — x')5{h — h'). 

The interaction potential Vr gives rise to a variety of phenomena. First, for ground-state 
problems the case of a attractive potential is obviously trivial: the two lines will localize 
to the same ground-state. In this paper we are going to deal with a hard-core interaction 
between the lines 1 and 2. This implies a delta-function-like Vr ~ Vo5{xi — X2)S{hi — 
with Vo — > oo so that overlap between the lines is strictly excluded. Would one allow for 
e.g. a finite Vq then the one-line ground-state would act as a pinning defect and the physics 
would slowly cross-over from the hard-core case to that of two independent lines as Vq is 
decreased. 

The simplest scaling picture for the TPRM in the presence of a hard-core interaction 
Vr consists of two independent directed polymers one being in the one-line global minimum 
and the second being in the first local minimum or the first excited state. This picture 
implies that the TPRM ground-state would be separable, that is it could be constructed by 
a successive optimization procedure. This turns out to be false, but the construction gives 
a definition for the effective interaction energy 

Vnt,eff = Ei+E[-E2r^L''' (2) 

where Ei refers to the single-DP ground-state energy in a particular sample, E[ to the first 
excited state, E2 is the true TPRM ground-state energy and L is the system size to be 



defined below in section Oy defines a scaling exponent for this particular form of the 
interaction energy. Recall that one has Ei ~ AL + AL^ + . . . and that the same is expected 
of E[ as well where A, A are disorder and dimension-dependent non-universal pre-factors. 
The argument is, however, essentially based on the claim that in the DPRM problem there 
is only one energy scale, that governed by the DPRM energy fluctuation exponent 6 and is 
therefore only qualitative. For E2 it is to be expected that the scaling is of the same form 
E2 ~ BL + BL^^ where the exponent 62 measures the energy fluctuations of the TPRM 



ground-state. The ensemble-averaged Vint,eff allows one to note that since the energy and 
its fluctuations have as an upper bound the separable trial ground-state 9v should be limited 
from above by 6i. 

Likewise, the interaction energy can be described by the energy of the TPRM ground- 
state minus twice the single line energy, i.e. 

5E2 = E2-2E1 L^^. (3) 

Here 9e defines another scaling exponent characterising the TRPM groundstate. One has 
naturally 5E2 + Vint,eff = E[ — Ei > and in particular if the single-line problem has two 
geometrically independent, energetically degenerate solutions then the sum is zero. Since 
6E2 is positive semi-definite sample-to-sample, a lower limit for 62 is 61 and therefore by 
this dual construction one would expect that 62 = Oi- In this work we do not consider the 
roughness properties of the two-line system but note that for it one would likewise expect 
that (2 = Ci- Figure shows examples from two and three dimensions of situations in 
which the true TPRM ground-state is non-separable, i.e. it can not be constructed out of 
the states with energies Ei and E[ and has thus a non-zero Vint,eff- 



III. NUMERICAL METHOD 

Here we define the lattice version of the continuum model of two random polymers with 
hard core interactions in a random environment introduced in the preceeding section. We 
formulate it in such a way that the connection to a minimum cost flow problem becomes 



obvious ||TT],|T2|, for which powerfull algorithms from combinatorial optimization exist that 
find exact ground states in polynomial time ||13|| . 
Consider the energy function 

H{x) = J2eij ■ x,j , (4) 

where is a sum over all bonds (ij) joining site i and j of a (i-dimensional lattice, 

e.g. a rectangular (L'^"^ x H) lattice, with periodic boundary conditions (b.c.) in d — 1 



space direction and free b.c. in one direction. The bond energies Cij > are quenched 
random variables that indicate how much energy it costs to put a segment of a polymer on 
a specific bond (ij). The variables describing the two polymers are Xij G {0,1} (for hard 
core interactions), Xij = 1 if there is a polymer passing bond {ij) and zero otherwise. For 
the configuration to form lines on each site of the lattice all incoming flow should balance 
the outgoing flow, i.e. the flow is divergence free 

V ■ X = , (5) 

where V- denotes the lattice divergence. Obviously the fiux-line has to enter, and to leave, 
the system somewhere. We attach all sites of one free boundary to an extra site (via 
energetically neutral arcs, e = 0), which we call the source s, and the other side to another 
extra site, the target, t as indicated in fig. la. Now one can push one line through the system 
by inferring that s has a source strength of +1 and that t has a sink strength of —1, i.e. 

(V ■ x), = +N and (V ■ x)^ = -N , (6) 

with = 1. Thus, the 1-line problem consists in minimizing the energy (^ by finding a 
flow X in the network (the lattice plus the two extra sites s and t) fulfilling the constraints 
(^) and (^. Naively one would expect that the 2- line problem consists simply in adding a 
second line to the 1-line configuration, avoiding the bonds already occupied due to the hard 
core interaction we consider here. A glance at Fig. |l| convinces us that this is not correct 
and actually the main issue of the present paper is to provide evidence that the correct 
solution of the TPRM problem is significantly different from what one gets when assuming 
the separability of the ground state. 

The first key ingredient to treat the two-line problem (and the A^-line problem in general 
pj| ) is that one does not work with the original network but with the residual network 
corresponding to the actual fiux-line configuration, which contains also the information 
about possibilities to send flow backwards (now with energy — Cjj since one wins energy 
by reducing Xij), i.e. to modify the actual flow. Suppose that we put one fiux-line along 
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a shortest path P{s,t) from s to t, which means that we set Xij = 1 for all arcs on the 
path P{s,t). Then the residual network is obtained by reversing all arcs and inverting all 
energies along this path, indicating that here we cannot put any further flow in the forward 
direction (since we assume hard-core interaction, i.e. Xij < 1), but can send flow backwards 
by reducing Xij on the forward arcs by one unit. This procedure is sketched in Figure 0. 

The second key ingredience is the introduction of a so called potential (p that fulfills the 
relation 

VU) < + eij (7) 

for all arcs (ij) in the residual network, indicating how much energy (p{j) it would at least 
take to send one unit of flow from s to site j, IF it would cost an energy ip{i) to send it to 
site i. With the help of these potentials one defines the reduced costs 

cf, = e,, + ^{{} - V9(j) > . (8) 

The last inequality, which follows from the properties of the potential (f (|^) actually ensures 
that there is no loop C in the current residual network (corresponding to a flow x) with 
negative total energy, since J2{ij)ec ^ij = J2{ij)(^c ^tji implying that the flow x is optimal []T2 



The idea of the successive shortest path algorithm is to start with an empty network, i.e. 
x° = 0, which is certainly an optimal flow for = 0, and set = 0, c^- = Cij. One now 
successively adds FL to the system using the following iteration: Suppose we have an optimal 
N — 1-line configuration corresponding to the flow x^^^. The current potential is v^^~^, the 
reduced costs are c^~^ = eij + Lp^~^{i) —'^^"^{j) and we consider the residual network G^~^ 
corresponding to the flow x^^^ with the reduced costs c^^^ > 0. The iteration leading to 
an optimal A^-line configuration is 

1. Determine shortest distances d{i) from s to all other nodes i with respect to the reduced 
costs cfj~^ in the residual network G^~^. 

2. For all nodes i update the potential: f^{i) = if^~^{i) + d{i) — d{t). 



3. Let P{s,t) denote a shortest path from node s to t. To obtain x^j increase (decrease) 
by one unit the flow variables xfj~^ on all forward (backward) arcs (ij) along P{s,t). 

(see Fig. ||). Note that due to the the fact that the numbers d{i) are shortest distances 
one has again > 0, i.e. the flow is indeed optimal. To estimate the complexity of this 
algorithm it is important to note that it is not necessary to determine shortest paths from s 
to all other nodes in the network; a shortest path from s to t is sufficient if one updates the 



potentials in a slightly different way ||T2||. Thus, the complexity of each iteration is the same 
as that of Dijkstra's algorithm for finding shortest paths in a network, which is (9(M^) for 
a naive implementation (M is the number of nodes in the network). We find, however, for 
the cases we consider (c?- dimensional lattices) it roughly scales linearly in M = L'^. Thus, 
for fiux-lines the complexity of this algorithm is 0{NL'^). 

In Figure § we show the true ground state configuration for a specific disorder configu- 
ration in 2d and in 3d and compare it with the one-line ground state plus the first excited 
state (the latter defined as the ground state in the network that is left when the bonds 
occupied by the one-line ground state are excluded). This is a typical example in which the 
two two-line configuration in 2d and in 3d are distinct. 



IV. RESULTS 

For the actual computations reported in the following we set the height of the system H 
equal to its lateral size L, i.e. H = L, yielding a square geometry in 2d and a cubic one in 
3d) and considered system sizes from L = 16 to L = 256 in 2d and from L = 8 to L = 64 
in 3d. For each system size the results are averaged over = 12000 (2D) and = 8000 
(3D) disorder configurations, and quantities like O = Ei, E2, SE2, Vint,eff denote disorder 
averages from now on. 

We expect the various exponents that we estimate to be independent of the actual disor- 
der we put in (as long it is uncorrelated and does not have algebraic tails), nevertheless we 
took two different probability distributions for the bond energies: 1) a uniform distribution 
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for which P{eij) = 1 for Cij G [0, 1] and otherwise; 2) a binary distribution in which Cij is 
1 with probabihty p and with probabihty 1 — p. 

A. Two dimensions 

Figure ^ shows the scahng of the two-hne system energy and energy fluctuations for 
both a uniform distribution for the Cj/s and a binary one with qp = 0.8. As expected, the 
scahng of the total energy E2 is linear and the fluctuations 6E2 scale with an exponent 62 
with $2 ^ 61, the one-line energy fluctuation exponent. This adheres to the picture that the 
energetics of the DP problem are in general dictated by the one-line exponent. 

In Figure ^ we show the probability that 6E2 = as a function of system size. This 
measures the true degeneracy of the two-line system as the joint ground-state can be obtained 
from two independent minima with the same energy. P{6E2 = 0) ~ L^"^ with ai = 
0.63 ± 0.03 which is compatible with to ai = 2/3 adhering thus to Tang's result |^ which 
indicated ai = 1 — 6. One can compare this with the scaling of P{Vint,eff = 0), which scales 
with an exponent 02 = 0.15 ± 0.02 for both distributions (P ~ L^"^). Similarly to Tang's 
conjecture, we are left with a picture which explains the frequency of separable ground-states 
(with 6E2 = 0) by a picture in which the two lines belong to two neighboring trees in the 
energy landscape. This means that one considers an inverted structure in which the two lines 
end up next to each other but belonging to two different trees (starting from x = L) with the 
same energy. Meanwhile the interaction energy in general shows increasing entanglement 
with a probability for a separable GS that decays with a novel exponent 02 = 0.15. 

Figures |^ and ^ discuss further the scaling of the mean interaction energies SE2 and 
Vint,eff for the both distributions. We find the exponents 9e ~ 0.39 ± 0.03 and 9v ~ 
0.39 ± 0.03, respectively. For both these quantities we seem to obtain that the effective 
scaling exponents are slightly higher than the one- or two-line energy fluctuation exponents 
as such. However, as shown in figure |^ we can collapse the energy probablity distributions 
for 6E2 and Vint,eff by using a two-exponent collapse. Note that this is different from the 
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simple collapse using 6e and By, however the two exponents combined make it so that the 
averages scale with 9e and 9v- 

B. Three dimensions 

Figure § shows the scaling of three-dimensional case again for both a uniform distribution 
for the Cij's and a binary one with p = 0.8 for the case of the two-line system energy 
and energy fluctuations. As expected, the scaling of the total energy E2 is linear and 
the fluctuations 6E2 scale with an exponent 62 with 62 — 61 ^ 0.24, the one-line energy 
fluctuation exponent in three dimensions. 

In Figure |^ we show the probability that 6E2 = as a function of system size. In 
3D, for binary disorder, P{SE2 = 0) ~ L^*^^ with ai = 0.25 ~ 61 in contrast with the 
geometric picture valid in 2D. The scaling of P{Vint,eff = 0) can not be described with a 
unique exponent and we find 02 = 0.11 ±0.01 for binary, and 02 = 0.05±0.01 for continuous 
disorder (P ~ L""^). Again the interaction energy in general shows increasing entanglement 
with a probability for a separable GS that decays with novel exponents ai, 02. 

Figures and |ll] discuss further the scaling of the mean interaction energies SE2 and 
Vint,eff for the continous distributions. The 3D exponents become 9e = 0.26 ± 0.02 and 
6v = 0.21 ± 0.02. As shown in figure |^ we can collapse the energy probablity distributions 
for 6E2 and Vint,eff by using as in 2D a two-exponent collapse. For binary disorder the 
collapse of the data makes sense in both cases, for continuous we restrict ourselves to 6E2. 

V. CONCLUSIONS 

In this paper we have investigated the joint ground-state of two directed polymers in 
a random medium, the TPRM problem. The main questions addressed here are whether 
the scaling of the TPRM can be described with the one-line exponents and an associated 
picture of behavior and if not so when. Unsurprisingly it turns out that Vint,eff as defined 
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here seems to result in an independent exponent that can not be explained by the one- 
line scaling arguments. This is natural since it measures the difference of the true TPRM 
ground-state to the 'Ansatz' of two separable states and is thus the first non-analytic and 
non-trivial correction characterizing the unique nature of the TPRM problem. On the other 
hand some of the features of the TPRM energetics, like the degeneracy of 5E2 are clearly 
related to the single-line picture in two dimensions. In three dimensions this is no longer 
true. We lack a geometrical explanation for the scaling of the degeneracy exponent 02 in 
this higher-dimensional case. 
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FIG. 1. Sketch of the successive shortest path algorithm for the solution of the minimum cost 
flow problem described in the text, (a) Network for = 0, the numbers are the reduced costs 
Cij, red for downward and right arcs, blue for upward and left arcs. The red line is a shortest path 
from s to t. (b) with the updated node potentials, (c) with the updated reduced costs. 
The green line is a shortest path, (d) OPtimal flow configuration for N = 2. 



13 



FIG. 2. a): Two-polymer ground-state in 2D, b) the same system but with the first (1-hne GS) 
frozen first, c) the TPRM GS in 3D, d) as b) but in 3D. In both the 2D and the 3D comparisons 
the disorder landscape is the same. 
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FIG. 3. Energy E2 (circles) and energy fluctuations AE2 (squares) of the TPRM problem in 
two dimensions in a log-log plot. We show data for binary disorder (cij £ {0, 1}) (filled symbols) 
and the uniform distribution of Cjj's (open symbols). One expects E2 (x L and AE2 oc , 
correspondingly the straight lines have slopes 1 (top) and 1/3 (bottom). 
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FIG. 4. P{5E2 = 0) (squares) and P{Vint,eff) = (circles) vs. L in 2d in a log-log plot. Data 
for both binary (filled symbols) and uniform (open symbols) distribution of the bond energies Cij. 
The data follow the relations P{6E2 = 0) oc L^"-^ and P{Vint,eff) = oc L""^ with ai and 02 given 
by the slopes of the straight lines: ai = 0.63 {^1 — 6) and 02 = 0.15. 
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FIG. 5. 6E2 in 2d for binary (filled symbols) and uniform (open symbols) distribution of the 
bond energies Cij in a log-log plot. It is 6E2 oc L^^ with 6e = 0.39 it 0.02. 
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FIG. 6. Vint,eff in 2d for binary (filled symbols) and uniform (open symbols) distribution of 
the bond energies e^j in a log-log plot. It is Vint,eff oc L^^ with 0^0. 39 ± 0.02. 
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FIG. 7. Scaling plots of the probability distributions of 6E2 and Vint,eff in 2d for binary (filled 
symbols) and uniform (open symbols) distribution of the bond energies Cjj in a log- log plot. 



16 



lO' 




FIG. 8. Energy E2 (circles) and energy fluctuations AE2 (squares) of the TPRM problem in 
three dimensions in a log-log plot. We show data for binary disorder (cjj G {0, 1}) (filled symbols) 
and the uniform distribution of e^j's (open symbols). One expects E2 oa L and A£^2 oc , 
correspondingly the straight lines have slopes 1 (top) and 0.24 (bottom). 
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FIG. 9. P{5E2 = 0) (squares) and P{Vint,eff) = (circles) vs. L in 3d in a log-log plot. Data 
for both binary (filled symbols) and uniform (open symbols) distribution of the bond energies eij. 
The data follow the relations P{SE2 = 0) oc L~"i and P{Vint,eff) = oc L~"2 with ai and 02 
given by the slopes of the straight lines: ai = 0.25 and 02 = 0.11, 0.05 for binary and continuous 
disorder, respectively. 
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FIG. 10. SE2 in 3d for binary (filled symbols) and uniform (open symbols) distribution of the 
bond energies Cij in a log-log plot. It is 6E2 oc L^'' with 9e = 0.26 ± 0.02 
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FIG. 11. Vint,eff in 3d for binary (filled symbols) and uniform (open symbols) distribution of 
the bond energies Cij in a log-log plot. It is Vint,eff with 9v = 0.21 ± 0.02. 
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FIG. 12. Scaling plots of the probability distributions of dE2 and Vint,eff in 3d for the uniform 
distribution of the bond energies e^j in a log-log plot. 
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